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The  Obukhov  dynamic  similarity  of  flows  hypothesis  for  the  surface  boundary 
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and  temperature  profiles  extracted  from  the  literature  led  to  the  conclusion 
that  no  constant  exists  for  the  log  plus  linear  profile  form  for  thermally 
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stratified  stable  flow  and  that  the  critical  value  of  the  gradient  Richardson 
number  is  unity. 
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INTRODUCTION 


Richardson  (1920)  originally  suggested  that  the  decay  of  turbulence  in 
the  atmosphere,  and  the  surface  boundary  layer  in  particular,  was  com- 
plete when  the  value  of  the  Richardson  number  Ri  became  unity.  This 
so-called  critical  Richardson  number  has  since  been  estimated  to  occur 
anywhere  in  the  range  of  0.0417  ■ Ri  <_  1.5.  A large  number  of  the 
estimates  of  the  critical  Richardson  number  have  been  based  upon  some 
form  of  a surface  boundary  layer  hypothesis  for  a thermally  stratified 
stable  regime,  with  Ri  ^ traditionally  inferred  from  the  constant,  . , 
of  the  ui aba  tic  influence  function.  There  have  been  nearly  as  many 
estimates  of  the  influence  function  constant  as  there  have  been  for 


Ri  . A cursory  review  of  vintage  and  contemporary  literature  re- 
veals that  0.6  • 17  generally  reflects  the  range  of  . , the  influence 

function  constant.  If  Ri^.^  = B_1,  then  the  values  for  Rifr^t  tabu- 
lated in  Table  1 as  determined  from  the  literature  review  are  in  good 
agreement  with  those  of  i-  given  in  Table  II. 

The  above  information  suggests  that  Ri ^ ^ should  be  near  unity 
and  that  is  more  likely  a variable  rather  than  a constant.  The  intent 


of  this  paper  is  to  explore  these  postulates. 

DYNAMIC  SIMILARITY  IN  A STABLE  REGIME 
Obukhov  (1946)  proposed  a dynamic  similarity  theory  for  the  surface 
boundary  layer  where  the  depth  of  the  surface  layer  was  characterized 
by  a unique  length  L defined  as 


L 


u*J^po 
kg  f~ 


(1) 
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TABLE  I 


Estimates  of  the  Value  of  the 
Critical  Richardson  Number  (0.667  _<  : _<  24) 


Source 

Critical  Ri 

Schlichtinq  (1036)  (See  Schlichtinq  1960) 

0.0417 

Rusinqer  (1955) 

0.105 

Sverdrup  (1936) 

0.11 

Obukhov  (1946) 

0.11 

McVehil  (1964) 

0.14 

Holzman  (194q) 

0.143 

Deacon  (194°) 

0.15 

Davis  (1957) 

0.15 

Ellison  (1957) 

0.15 

Lumley  and  Panofsky  (1964) 

0.14  to  0.22 

Busch  (1965) 

0.20 

Businger  (1959) 

0.20 

Webb  (1970) 

0.20 

Businger,  et  al . (1971 ) 

0.21 

f-tancon  Pi  p H ( ] OflQ  \ 

9.224 

Goldstein  (1931) 

0.25 

Ludlam  (1967) 

0.25 

Miles  (1961,  1963) 

0.25 

Sieffert  (1958) 

0.25 

Taylor  (1931 ) 

0.25 

Kaimal  and  Izumi  (1965) 

0.25 

Webster  (1964) 

0.25  to  0.50 

Businger  (1969) 

0.25  to  1.00 

Townsend  (1957) 

0.30 

Oke  (1970) 

0.30 

Wanta  (1953) 

0.33  to  0.50 

Portman  et  al.  (1962) 

0.35 

Crawford  (1965) 

0.35 

Rider  and  Robinson  (1951) 

0.40 

Lyons,  et  al . (1964) 

0.50 

Sutton  (1949) 

0.50  to  1 .00 

Fichtl  (1970) 

0.55 

Cal der  (1949) 

1 .0 

Petterssen  and  Swinbank  (1947) 

1 .0 

Durst  (1933) 

1.0 

Endlich  (1964) 

1 .0 

Kronebach  (1964) 

1 .0 

Richardson  (1920) 

1 .0 

Zavarina  and  Yudin  (1960) 

1 .0 

Clodman  (1953) 

1.5 

Jaffe  (1963) 

1 .5 

i 


TABLE  II 


Estimates  of  the  Linear  Influence 
Function  Constant  (0.058  £Ricr^t  ± 1.67) 


Source 
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Arya  (1968) 

Busch  (1965) 

Businqer,  et  al . (1971 ) 

Chalikov  (19687 
Deacon  (1955) 

Deardorff  ( i%8) 

Fichtl  and  McVehil  (1970) 

Gurvich  (1965) 

Hansen  and  Serna  (1969) 

Holzman  (1943) 

Lumley  and  Panofsky  (1964) 

McVehil  (1964) 

Monin  and  Obukhov  (1954) 

Oke  (1970) 

Panofsky  (1967) 

Panofsky,  et  al . (1960) 

Plate  and  Lin  (1966) 

Rossby  and  Montgomery  (1935) 
Sverdrup  (1936) 

Webb  (1970) 

Zil itinkevich  and  Chalikov  (1968) 


10  17 

5 

4.7  4.7 

3.77  to  7.44  5.1 

7.1  to  12 

7 11 

7 

8.5 

7 10 

7 
7 
7 

0.6  0.6 


5 

10 

4.5 

7 

11 

11 

5.2 

9.91  10.4 


*"  values  determined  from  temperature  profile  analysis 


4 


7 to  7.19 


where  u*  is  a friction  velocity,  c^  the  specific  heat  of  air  at  constant 
pressure,  . density,  e potential  temperature,  k Karman's  constant,  g the 
yra vi tational  acceleration,  and  H the  vertical  heat  flux.  The  scaling 
length  L evolved  from  the  integration  of  the  Richardson  number  with  respect 
to  height  as  a function  of  the  dimensionless  wind  shear  and  lapse  rate. 

From  the  Obukhov  hypothesis,  the  deduction  may  be  made  that  the  wind  and 
temperature  profiles  in  differential  form  are  given  by 


aV 

jZ 


U* 

kz 


(2) 


and 


where  V is  the  mean  horizontal  windspeed,  z is  height,  0,^ 
dimensionless  shear  and  lapse  rates,  respectively,  and  T* 
temperature  given  by 


(3) 

and  0,,  are  the 
is  a scaling 


T* 


11  _ 
CP  P 


(4) 


The  eddy  viscosity  and  eddy  conductivity  are  then  generally  written  as 


km  = ku*z  0;V 


(5) 


and 

Klt  = ku*z  r\\  ■ (6) 

If,  in  a stable  regime,  the  transfer  of  heat  and  momentum  are 
assumed  to  be  a function  of  mechanical  turbulence  only,  as  may  be  in- 
ferred from  Reynolds  analogy  and  the  vertical  turbulent  kinetic  energy 
budget,  generally  stated  as 


I|2  '±1  + flH  . i§w 
* jz  Cp  pe  jz 


iwp/p0 

JZ 


- t = 0 


(7) 
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T T 


where  e is  total  kinetic  energy,  p pressure,  w the  mean  vertical  component, 
and  the  viscous  dissipation  rate,  then 


and  by  defi ni tion 


KH  KM;  0M  ' 


f = Ri  0M 


(8) 


where  z/L  is  the  familiar  Mon i n and  Obukhov  (1954)  scaling  ratio.  From 
Equation  (8),  L may  be  redefined  as 


u*G3V/az 

L * 


kgae/  )z 


R-j  = g as/  az  i 
e (aV/3z) 


(9) 


(10) 


with  0,_,  defined  by  Equation  (2). 
Rewriting  Equation  (2)  as 


aV  . u*  z/L 


(11) 


az  kz  Ri 

then  adding  and  subtracting  1 to  z/L  (Ri)-1  and  multiplying  and  dividing 
by  z/L  yields 


IV  = * i + L h/k  - Ri 
32  kz  L I 81  z/L 


(12) 


Defining  (z/L  - Ri ) (Ri  z/L)  ' as  an  arbitrary  variable  i- , then  inte- 
grating Equation  (12)  gives 


t) 


(13) 


where  z^  is  the  roughness  length  and  IT  is  the  average  . over  a layer 

Z = 


I 


It  is  easily  demonstrated  that  z/L  0r_^  - 1 so  that  Lquation  (13)  may 
be  stated  as 


" f-  * <D„  - 1)  . 

0 


If  0.^  = 0((,  then  the  inteyrated  temperature  profile  is 

6 - en  = T*  in  f-  + (0M  - 1)  (lb) 

u L 0 

where  ;,-g  is  the  potential  temperature  at  z^. 

Evaluation  of  the  profile  hypothesis  is  dependent  upon  establishing 
L from  experimental  data.  Equation  (9)  shows  that  the  only  inferred  entity 
is  the  friction  velocity  u*.  All  other  parameters  are  constant  or  can  be 
easily  measured.  The  friction  velocity  can  be  evaluated  by  first  consid- 
ering the  inteyrated  form  of  the  wind  profile  [Lquation  (13)].  If  the 
mean  wind  \l  l at  Z[  is  subtracted  from  speed  V 3 at  z3,  where  the  subscripts 
1,  2,  3,  ...  refer  to  adjacent  levels. 


u*  Z3  _ z3 

-1  = \T  77+  e3  r 


- Z1  I 

■ r r 


In  a geometric  progression  such  as  y = ar  with  a = 1 and  r = 2,  the  rela- 
tionships between  the  levels  may  be  expressed  as  z3  = 2z  and  zl  = 1/2  z . 


Thus,  Lquation  (16)  can  be  rewritten  as 


'3,1  “ 


- 

z.i  / - 6l  \ 

+ r y ^ /'  • 


From  the  profile  di fferential  , the  following  can  be  determined: 


z3  - z: 


,*  = F ^ - + ^n  - e2  ri  • 
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From  Equations  (17)  and  (18) 


4i- 


~r 


3 _ 

'n  7,  ez‘ 


(19) 


Substitution  of  Equation  (19)  in  (17),  after  recalling  that  V = u*/k 
[■n  z2/zQ  + (3 2 z2/L]  and  that  for  a given  profile  u*/k  is  constant, 
gives 


a2 

r 


tu 


3,1 


*2  z 

in  r V2  in  - 


z2  V2  an 


L0 
Z 3 

:1 


(20) 


The  friction  velocity  for  each  profile  may  be  found  from  Equation  (13), 
if  the  roughness  length  Zq  is  known.  For  the  two  data  samples  used  in  the 
evaluation,  values  of  z^  were  based  upon  those  determined  by  Panofsky  (1963). 

VALIDATION  OF  THE  HYPOTHESIS 


The  hypothesis  was  validated  by  using  thermally  stratified  stable  regime 
data  extracted  from  the  Great  Plains  Turbulence  Program  (Lettau  and  Davidson 
1957  and  Project  Prairie  Grass  (Barad  1958).  All  stable  profiles  (Ri  > 0) 
were  plotted  as  V,  e,  = f (nn  z)  and  inspected  for  suspect  erroneous  data 
points.  The  vertical  gradients  were  determined  such  that  the  gradients 
were  tangent  to  the  geometric  mean  heights  of  the  layers.  Of  the  163 
profiles  available,  103  were  retained  for  use.  Friction  velocities  and  the 
scaling  length  L were  calculated  for  each  profile.  Fluctuations  in  the 
values  of  L with  respect  to  height  were  smoothed  by  using  an  approximation 
of  the  height  derivative  in  the  form 


8 


r 


Values  of  the  dimensionless  shear  were  calculated  from  equation  (2). 

The  Richardson  numbers  for  the  profiles,  averaged  over  suitable  in- 
tervals, are  shown  as  a function  of  z/L  in  Figure  1.  The  bars  represent 
one  standard  deviation  about  the  means.  The  general  characteristics  of 
the  profile  shape  indicated  that  z/L  as  a function  of  Ri  was  parabolic. 
Consequently,  the  experimental  data  were  fitted  by  least  squares  methods 
to 

^ = a + b Ri  + c Ri'  (22) 

which  leads  to  a = 0,  b = 1,  and  c = 15,  so  that 

f-  = Ri  + 15  Ri?  (23) 

The  solid  curve  in  Figure  1 represents  Equation  (23). 

Values  of  .•  for  each  profile  layer  were  calculated  from 


. z/L  - Ri 


(24) 


Ri  z/L 

After  averaging,  . as  a function  of  z/L  and  Ri  was  plotted  as  shown  in 
Figure  2.  From  Equations  (7)  and  (23), 


7 - 15 

3 " TT  TeTrI 


(25) 


and 


0H  = 1 + 15  Ri 


(26) 


The  curve  in  Figure  2 represents  Lquation  (25). 
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THE  CRITICAL  RICHARDSON  NUMBER 


Richardson's  (1920)  criterion,  in  terms  of  the  production  of  energy  can 
be  shown  to  be 


dE 


U,t) 

~3t 


/ — \2 

1 3V 

Si  g 

ae/ez 

\ 32  / 

> " 6 

bV/oz)2 

(27) 


where  dE(z,t)/dt  is  the  rate  of  increase  of  the  mean  turbulent  kinetic 
energy  per  unit  time.  Factoring  K^/ and  0M  from  the  right-hand  side 
of  Equation  (27)  yields 


dEu.t)  . Si  M 

dt  0M  \zz  ) 


(28) 


which  implies  that  the  complete  decay  of  turbulence  does  not  occur  until 
z/L  approaches  the  value  of  the  concurrent  0^.  From  Equations  (23),  (25), 
and  (26),  critical  0M  and  z/L  appear  to  be  about  16,  which  can  also  be 
inferred  from  Figures  1 and  2.  This  strongly  suggests  that  Ri  ^ is 
unity.  Values  of  R1 , 0^,  and  e as  a function  of  z/L  as  calculated  with 
Equations  (23),  (25),  and  (26)  are  given  in  Table  III. 


DISCUSSION  AND  SUMMARY 

The  characteristics  of  thermally  st>;tified  stable  flow  as  deduced  from 
the  Obukhov  (1946)  similarity  theory  and  inferred  from  the  vertical  tur- 
bulent kinetic  energy  budget  should  account  for  the  wide  range  of  critical 
values  of  the  gradient  Richardson  number  reported  in  the  literature. 

Many  values  of  R1  ^ were  evidently  based  upon  values  of  e'1  calculated 
from  data  observed  over  an  extremely  narrow  range  of  Ri  or  z/L.  The 
values  of  e extracted  from  the  literature  as  a function  of  z/L  are  shown 
in  Figure  3.  Again,  the  implication  that  u is  variable  is  quite  apparent. 
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TABLE  III 


The  Richardson  Number,  Dimensionless  Shear,  and 
. as  a Function  of  Z/i  for  Stable  Flow. 


z/i 

Ri 

8 

0 

0 

1 

15 

0.0012 

0.001 

1 .0150 

14.78 

0.0025 

0.002 

1 .0300 

14.56 

0.0054 

0.005 

1 .0750 

13.95 

0.0077 

0.007 

1 .1050 

1 3.57 

0.0115 

0.010 

1 .1500 

13.04 

0.0260 

0.020 

1.3000 

11.54 

0.0435 

0.030 

1 .4500 

10.34 

0.0640 

0.040 

1 .6000 

9.38 

0.0875 

0.050 

1.7500 

8.57 

0.10 

0.055 

1 .8250 

8.22 

0.20 

0.087 

2.3050 

6.51 

0.30 

0.1120 

2.6800 

5.60 

0.40 

0.1335 

3.0030 

5.00 

0.90 

0.1525 

3.2875 

4.56 

0.60 

0.1695 

3.5425 

4.23 

0.70 

0.1860 

3.7900 

3.96 

0.80 

0.2000 

4.0000 

3.75 

0.50 

0.2140 

4.2100 

3.56 

1.0 

0.2270 

4.4050 

3.41 

2.0 

0.3340 

6.0100 

2.50 

3.0 

0.4154 

7.2310 

2.07 

4.0 

0.4845 

8.2675 

1 .81 

5.0 

0.5450 

9.1750 

1 .63 

6.0 

0.6000 

10.0000 

1 .50 

7.0 

0.6507 

10.7605 

1.39 

8.0 

0.6976 

11.4640 

1 .31 

9.0 

0.7420 

12.1300 

1.29 

10.0 

0.7840 

12.7600 

1 .18 

11.0 

0.8290 

13.3600 

1.12 

12.0 

0.8620 

13.9300 

1 .08 

13.0 

0.8985 

14.4775 

1.04 

14.0 

0.9335 

15.0025 

1 .00 

15.0 

0.9675 

15.5125 

0.97 

16.0 

1 .0000 

1 6 . 0000 

0.94 
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FUNCTION  OF  z/L  AS  EXTRACTED  FROM  THE  LITERATURE 


The  concepts  explored  in  this  treatise  are  all  straightforward  and 


generally  agree  with  the  mean  flow  conditions  observed  in  the  surface 
boundary  layer  in  thermally  stratified  stable  air.  The  results  are  also 
in  reasonably  good  agreement  with  the  results  of  other  experimenters,  if 
the  linear  influence  function  and  Ri  of  these  studies  are  adjusted 
to  fit  current  results.  The  vast  range  of  critical  gradient 
Richardson  numbers  and  the  "constant"  £ found  in  the  literature  tend 
to  support  the  author's  contentions.  An  interesting  sidelight  of  this 
study  is  the  adiabatic  value  of  £ = 15,  which  coincides  with  the  value 
of  y = 1 5 for  the  linear  quartic  form  of  the  unstable  profile,  usually 
written  as  0^  = (1  - y z/L)"‘/‘,  as  reported  by  many  experimenters. 

CONCLUSIONS 

Two  conclusions  may  be  stated:  (1)  The  Critical  Richardson  number  is 
unity,  and  (2)  no  constant  exists  for  the  log-plus-linear  influence 
functi on. 
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